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The polarization operator of an axion in a degenerate gas of electrons 
occupying the ground-state Landau level in a superstrong magnetic field H ^ 
' Hq = m^c^/eh = 4.41 • 10^^ G is investigated in a model with a tree-level 

a^ : 

' axion-electron coupling. It is shown that a dynamic axion mass, which can 

fall within the allowed range of values (10~^ eV < nia ^ 10^^ eV), is generated 
\^ ' under the conditions of strongly magnetized neutron stars. As a result, the 

dispersion relation for axions is appreciably different from that in a vacuum. 

> : 

^ I 1. The a priori strong nonconservation oi CP parity in the standard model can be elim- 

^ I inated in a natural manner by introducing axions — pseudo-Goldstone bosons associated 

^ I with the spontaneous breaking of the additional Peccei-Quinn global symmetry U{1)pq. 

QhI According to the experimental data,i the energy scale Va for the U{1)pq symmetry breaking 
Oh' 

Q I is much greater than the electroweak scale — 10 GeV, and the constants of the pos- 



sible couplings of an axion to the standard particles (~ 1/fa) are very small (the "invisible" 



^ I axion: see Ref. | for a review of various axion models). 



Axion effects can be appreciable under the astrophysical conditions of high matter den- 
sities, high temperatures, and strong magnetic fields (for example, in neutron stars!) . Axion 
production processes, which result in additional energy losses by stars, and the limits ob- 
tained by astrophysical methods on the parameters of axion models are examined in Ref. ^. 
In so doing, the influence of electromagnetic fields was neglected. 

The investigation of axion processes in strong magnetic fields commenced comparatively 
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recently. The Compton and Primakoff mechanisms of axion production on nonrelativistic 
electrons by thermal photons (7 + e — > e + a) in the presence of a magnetic field are studied 
in Ref. ^. The extension to relativistic electrons in a constant external electromagnetic field 
is given in Ref. ^ (Primakoff effect) and Refs. || and ^ (Compton effect), whereBi estimates 
were also obtained for the contributions of the indicated processes to the axion luminosity of 
a magnetized strongly degenerate relativistic electron gas under the conditions of the crust 
of a neutron star. A new axion production mechanism — synchrotron emission of axions 
(e — s> e + a) by relativistic electrons — was proposed in Ref. |10| and its contribution to the 
energy losses by a neutron star was calculated. In Refs. ^ - |l^ it was assumed that the 
external field intensity F <^ Hq = rn^c' /eh ~ 4.41 ■ 10^^ G. In Ref. |ll], numerical methods 



were used to extend the results of Ref. to superstrong magnetic fields H > Hq. It was 
found that the basic equation derived in Ref. |1^ for the axion synchrotron luminosity for 
the semiclassical case of high electron energies (e ^ mgC^) and fields H <^ Hq agrees with 
the numerical calculations up to H/Hq < 20. The axion synchrotron luminosity of neutron 
stars and white dwarfs was also investigated in Ref. |TT[ 

In Refs. ^ - [IT] a model with a derivative axion-electron coupling eae, described by the 
interaction Lagrangiani 

^ae = 1^ (V^7^V) d,a, (1) 

was used. Here rrie is the electron mass and 7^ = —i'y'^'-f^'-f'^'y^; the system of units such that 
h = c = 1 is used; the signature of the metric is (H ); and 

Qae = Ce (2) 

is a dimensionless coupling constant, where the numerical factor Ce depends on the choice 
of the specific axion model.S 

In models where axions are coupled only with heavy fermions by a tree- level coupling 
there arises an effective direct low-energy axion-photon interaction of the type 707 .@ This 
interaction is the basis of the Primakoff axion photoproduction mechanism employed in 



Refs. ^ and ^ The synchrotron process e — ea in the absence of a tree-level axion-electron 
coupling was considered recently in Ref. [1^. This process is due to resonant conversion 
of a longitudinal plasmon (a photon in a medium), emitted by a relativistic electron in a 
magnetic field, into an axion. 

Decay of an axion in a strong magnetic field into a fermion pair (a — > and two 

photons (a 77)^ are also of interest for astrophysics and cosmology. 

In the present paper the model (1) is used to calculate the polarization operator of 
an axion moving in a strongly magnetized degenerate electron gas and the change in the 
dispersion relation of an axion in a medium is investigated using this operator. 

2. Taking account of the contribution of the electrons only (see Eq. (1)) we obtain, using 
the real-time formalism of the finite-temperature quantum field theory (see, for example. 



Ref. |T5D, the following momentum representation for the one- loop polarization operator of 



an axion: 



Il{k,k') = —iG^ J d X d x' exp{ikx — ik'x')TT kj G{x,x')k''y G{x' ,x) . (3) 

Here k {k') is the final (initial) 4-momentum of an axion; G{x,x') is the time-dependent 
single-particle Green's function of an ideal electron-positron gas in a constant magnetic 



field;El notations have also been introduced for the contraction a = 7'^a„ of a 4-vector 
with the Dirac 7 matrices and for the dimensional coupling constant 

On account of the translational invariance (constant external field, homogeneous isotropic 
medium) the polarization operator (3) is diagonal in momentum space: 

n(fc, k') = (27r)V^)(fc - k')U{k). (5) 

Here H(fc) determines the axion propagator D{k) in the momentum representation according 
to the Dyson equation 

Dik)=[e-ml-Uik)]-\ (6) 
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where is the free-axion mass (in the absence of a field and a medium), which is generated 
by the chiral anomaly of QCD:I rria ~ h?Q(jj^/va- The renormalized value njj(/c) (see below) 
gives the dispersion relation 

k^=ml + YiR{k). (7) 

3. We give the constant uniform magnetic field H || z in terms of the 4-potential in 

the gauge 

A'^ = (0,0,a;/7,0). (8) 

Then the Green's function G{x,x') can be represented in the following form after summing 
over the spin quantum number and the sign of the energy in the general expression for G in 
the form of a series in quadratic combinations of the eigenfunctions of the Dirac operatorlll: 

G{x, x') = [Y{id^ + eA^) + rrie] K{x, x'); 

oo 

K{x, x') = — — ^ ^ / dpodpydp^ exp [-ipo{t - t') + ipy{y - y') + ip^{z - z')] x 
^^^'^ „=o 



-oo 



X UnivWW) {Rn+l^+ + Rn^-) , (9) 

Rn = [pI -pI- '2hn -ml + iO] + 2ni6{pl -pi- 2hn - ml)NF{po). 

Here the electron charge —e<0,h = eH; n = 0, 1, 2, ... is the principal quantum 
number (the number of the Landau level); Py and Pz are the eigenvalues of the operators of 
the canonical momenta — the constants of motion in the gauge (8); and Unirj) is a Hermite 
function of argument 

r] = Vh{x+Py/hy, 7]' = r]{x x'); S± = (1 ± S3) /2, S3 = ^tV- 

The first term in i?„ has poles at the points po = = ±[m1 + 2hn+p1Y^'^, determining 
the energy spectrum of an electron in a magnetic field. The second term (oc 6{pq — e"^)) 
describes the effect of the electron-positron medium, and 



Nf{po) = 0{po) [exp [P{po - fj,)] + 1]-' + 9{-po) [exp [P{-po + fx)] + 1]'^ (10) 

is expressed in terms of the Fermi distribution function of electrons and positrons in a 
medium with temperature T — 1/(3 and chemical potential /i, and 9{±po) is the Heaviside 
step function. 

4. It is difficult to make a general analysis of the axion polarization operator for arbitrary 
values of the parameters H, T, and /i. In the present paper we confine our attention to 
superstrong magnetic fields and comparatively low temperatures 

H:^Ho, T<^n-me, (11) 

and we require the chemical potential to satisfy 

11^ -ml< 2h. (12) 

It follows from Eqs. (11) and (12) that in this case the contribution of positrons in Eq. 
(10) can be neglected (it is suppressed by the factor exp[— /?(// — me)]) and the medium is a 
degenerate gas of electrons occupying the ground-state Landau level (n — O): 

NF(po)^e(po)e(ii-po), Po^Vml+pl (13) 

We also limit the range of the axion 4-momentum 

\kl-kl\<^h. (14) 

Then the main contribution of virtual (vacuum) electrons and positrons is likewise formed 
by states with n = 0. As a result, retaining on the basis of Eqs. (11), (12), and (14) terms 
with n = in the sum (9), we obtain the following approximate expression for the Green's 
function in a superstrong magnetic field: 



— oo 

X / [-'Poit - t') + 'P^i^ - ^')] G{p)E.. (15) 



-^{v^ + v'^) + ipyiy-y') 



Here p = {po, 0, 0,pz) and 



G{p) = {p + me) [(p^ -ml + iO)-^ + 2-ni5(j)'^ - ml)NF{po)\ 



(16) 



is the Fourier transform of the Green's function in the two-dimensional space (0, 3). For 
Np = (no medium) the expression (15) is the well-known, effectively two-dimensional, 
electron propagator used in the theory of electrodynamic processes in superstrong magnetic 
fields and, specifically, for investigation of the photon polarization operator. 

5. Let us substitute the expression (15) into Eq. (3) and integrate over t,t' ,y,y', z, and 
z'. This gives in the form of a product of delta functions 

Y[ ^{K- kn)S{Pn + kn- Pn) 
n=0,y,z 

the laws of conservation of energy and of the corresponding projections of the momentum. 
The subsequent calculation of the Gaussian integrals over x, x' and the trivial integral over 
Py gives 6{k'^ — k^). As a result, as should be the case, we obtain a diagonal representation 
of the polarization operator (5), where 



U{k) = ^hexp (-^][F{l) + M{l)] 



(17) 



F{1) = 



(18) 



M{1) = 27r 



(Pp 



Sip^ - mDNpipo) 



ni,p) 



{p — ly — ml + iO 



+ -I) 



Here p = {po, 0, O,^^) and / = (/cq, 0, 0, kz) are two-dimensional vectors, and 
T{l,p) = ^Tr A;7^(p + me)S_^7^ (jj - I + m^^ S_ 



(19) 



(20) 



In Eq. (17) the function F corresponds to the purely field contribution, and M describes 
the influence of the medium. We note that M does not contain a term oc Np{pq)Np{pq — Iq), 
since 
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Sip' - ml)5 {(p - If - ml) 0ipo)eipo - k) = 0. 

Using the relations 

[E_,p] = 0, [S_,7^] = 0, 7"S_ = E+7" (n = l,2), E+E_ = 

the trace of Eq. (20) reduces to a two-dimensional form and can be easily calculated as 

T{l,p)^^Tr [f(p + me)f(p-Z-me)] ^ 2{lpf - f{lp + + ml). (21) 

We calculate the Gaussian integrals over po and Pz in Eq. (17) using the trace (21) and 
the well-known Fock-Schwinger proper-time representation for propagators of the form 



(A + iO)"^ = -i J ds exp [is{A + iO)] . 





As a result, we find for the function F{1) = Fip^^m^) the integral representation 

^1 oc 

F{P) = -i^ jdv j dx {[1+{1-v')t] exp [-ix [1-(1-^;')t]] -exp(-ix)} , (22) 



T = f/Aml. 

Here regularization is performed according to the well-known rule 

F{f) = F{f,m')-F{f,K') 

with the regulator mass A — > oo. 

The renormalized polarization operator ^r{t) can be obtained ^from the dispersion 
relation with one subtraction (as well as it was done for the photon polarization operator in 
Ref. 16) 

oo 

in,(r) = i / f"-"""' (23) 

T ^ ' J t{t-T -iQ) ^ ' 



Eq. (22) gives immediately 
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ImF = y dv5 v^)t] = -mlQ{T -\)[\ - 



-1/2 

(24) 



For r < we obtain from Eqs. (23), (24), and (17) the field contribution 

n<f..-^,„J^exp(4)i^K (25) 

to the axion polarization operator. Here = gl^./'^'n (see Eq. (4)), and the standard 
varia ^ was introduced as 

which is convenient for analytical continuation in = AitlIt. 

For r > 1, a channel is open for axion decay into an electron-positron pair [a —>■ e^e^) 
in a magnetic field. Its rate w for a real axion is related with the imaginary part of the 
polarization operator on the mass shell by the well-known relation 

.^-il™n<f'.a„!^^e.p(-|).(.-l)(l-i)-", (27, 

where w is the axion energy. 

This result, which follows from Eq. (25) with ^ = |^| exp(27r) (see Eq. (26)), is identical to 
the result obtained in Ref. ^ on the basis of a calculation of the elastic scattering amplitude 
of an axion in a magnetic field. It can also be found at once from Eqs. (24) and (17). 

Let us consider the contribution M (19) of the medium to the axion polarization operator. 
We note that it does not renormalize.El Integrating over in Eq- (19); using the delta 
function and taking account of Eqs. (13) and (21), gives 

M = ^ [D{l,p) + D{-l,p) + Dil,p) + Di~l,p)] , 

In J q 

me 

D{l,p)=[l^-2{lp) + i0]~\ (28) 



Here e is the energy of electrons in the medium, q = a/^^ — rnl, and the two-dimensional 
scalar products are Ip = k^e — k^q and Ip = k^e + k^q. 
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The imaginary part of the expression (28) is determined using Sokhotskii's formula 

^ -^P- -i7r5{x), (29) 



X + iO X 

where P signifies a principal value. From Eqs. (28) and (29) we obtain on the mass shell 

2 

ini M = -fe{T - 1) [^(// - £+) + eifx - £_)] , 

Here e± are the roots of the equations — 2uje ± 2kzq = 0. 
^From Eqs. (17), (24), and (30) we find the rate 

= ^ [0{e+ - IJ) + 9{e- - n)] w, (31) 

where w is the decay rate (27) in the absence of a medium, for the axion decay into an 
e~e'^ pair in the presence of a magnetized degenerate electron gas. We underscore that the 
imaginary part of the contribution (30) of the medium is positive, and summed with the 
negative field contribution (24) it gives a blocking Pauli factor 1 — 9{x) = 9{—x) in Eq. (30). 
It forbids electron production inside a filled Fermi sphere (for e± < ii). 

Taking account of Eq. (29), we obtain for the real part of Eq. (28) on the mass shell the 
representation 



A 

771 ^ 

ReM = ^-r/dx 

TT 



+ 



(32) 







_T — cosh^(x — ■0) T — cosh^(x + '0)_ 

Here the substitution of the variable £ — > x was used: e — mg cosh x and q — nie sinh x, and 
the parameters A and ip, defined as 

cosh A = — , tanh'^ = — , (33) 
rrie cu 

were introduced. The integral (32) can be expressed in terms of elementary functions. 

We shall confine our attention below to the limiting cases that are of interest for astro- 
physical applications. 



6. For an axion on the mass shell 



/2 = 4mV = 



ml + k] > 0, 



(34) 



and the condition (14) gives k\ <C h, so that exp(— /c^/2/i) ^ 1. We note that the imag- 
inary part of the polarization operator is formed by the contribution of real electrons and 
positrons, and the expression for it holds under the weaker condition k'j_ < 2h. Therefore 
the exponential factor can be retained in Eq. (27). 

For r ^ 1 (substantially below the threshold of the decay process a e~e~^), we find 
from Eqs. (25), (32), and (17) 



n 



R 



H 



(35) 



Here 



v± = tanh(A ± ip) 



V = tanh A 



1-f^ 



1/2 



Lj ± ukz 

We note that < and if the axion moves in the direction of the field H {kj_ = 0), then 
according to Eqs. (34) and (35) Ur — > in the limit of a massless axion {rria 0). 

At high energies (r ^ {fi/rrif,)'^ ^ 1) we obtain for the polarization operator the asymp- 
totic representation 



TT Ho 



In 



2aa 



eH 



IT 



In t— 

2fi 2 



(36) 



and it does not depend of the electron mass rrig as it should be in this limit. 
Let us write the dispersion relation (7) in the form 



kl + kl + ml + nR{k). 



(37) 



It follows from Eqs. (35) - (37) that in a magnetized medium a radiative shift of the axion 



mass is generated — a dynamic mass, whose square, according to the definition in Ref. |T^, 
is 



Smi = Re n 
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For r > 1 and /i/mg ^ 1 we obtain the estimate 

5ma gaemj—j r^lO^g^l-^j eV. (38) 

For gae ~ 10-13 (^j^gfg | ^^^j |10|) and i7 > 10^^ G (such fieldJlllli and even H ~ 10^^ - 10^° 
G (Ref. |2y) can exist in the interior regions of neutron stars), Eq. (38) gives Snia > 10"^ 
eV. 

The chemical potential /i of a degenerate gas of electrons occupying the ground-state 
Landau level (n = 0) in a magnetic field is related with the electron density Ue by the 
well-known relation 

where Pf = \/ — "ml is the Fermi momentum. Writing Eq. (12) in the form 



2 



->-f-), (40) 
Ho 2\mJ ' ^ ^ 



we obtain, taking account of Eq. (39), an upper limit on the density 

(41) 



where Xe = 1/me is the electron Compton wavelength. For H = 2-10^^ G Eqs. (40) and (41) 
give Pf < 50 MeV and < 10^^ cm'^. Next, let T ~ 10^° K ~ 1 MeV and > T. Then 
the conditions (11), (12), and (14) can be satisfied and the estimate (38) can be justified. 

In summary, under the conditions of strongly magnetized neutron stars a dynamic axion 
mass, which can fall within the existing limits on the axion massi'SIll — 10~^ eV < rria < 10^^ 
eV — is generated. Therefore Srria ~ rria and the dispersion relation (37) differs appreciably 
from the vacuum relation (fc^ = m^). This must be taken into account, for example, when 
investigating the resonant conversion of a plasmon into an axion (7 — > a) in a magnetic field 
as a result of the crossing of the corresponding dispersion curves (as already noted above, 
this process in fields H <^ Hq and in the absence of the direct coupling (1) was studied in 



Ref. |T2|). We also note that the rate (27) of the decay a e e"*" in a magnetic field has 
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a square-root threshold singularity (as r — > 1 + 0). This singularity can be removed by 
taking into account accurately the dispersion law of an axion near threshold, and the decay 
rate is found to be finite:lii w ~ me{ciaH / Hq)"^/^ . A detailed analysis of the same threshold 
singularity (of cyclotron resonance) in a magnetic field and its elimination for the photon 
decay process (7 — » e^e+) was given earlier in Ref. where, specifically, it is underscored 
that the indicated singularity can be explained by the quantization of the phase space of 
charged particles in a magnetic field. 
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